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We investigate the structure and features of an ultralong-range triatomic Rydberg molecule formed 
by a Rb Rydberg atom and a KRb diatomic molecule. In our numerical description, we perform a 
realistic treatment of the internal rotational motion of the diatomic molecule, and take into account 
the Rb(n, I > 3) Rydberg degenerate manifold and the energetically closest neighboring levels with 
principal quantum numbers n' > n and orbital quantum number I < 2. We focus here on the 
adiabatic electronic potentials evolving from the Rb(n, I > 3) and Rb(n = 26 ,1 = 2) manifolds. 

The directional properties of the KRb diatomic molecule within the Rb-KRb triatomic Rydberg 
molecule are also analyzed in detail. 


I. INTRODUCTION 

Recently the existence of a new class of ultralong-range 
polyatomic Rydberg molecules was predicted theoreti¬ 
cally [3, i- These huge-sized Rydberg molecules are 
formed by a Rydberg atom and either a A-doublet or a ro¬ 
tating polar diatomic molecule EIJ. The binding mech¬ 
anism is established by the electric field of the Rydberg 
electron and the ionic core that couples either the two 
internal states or hybridizes the rotational motion of the 
molecule due to its permanent electric dipole moment. 
The existence of these ultralong-range triatomic Rydberg 
molecules was predicted for polar diatomic molecules 
with subcritical electric dipole moments (do < 1-639 D) 
in order to prevent the binding of the Rydberg electron 
to the heteronuclear (polar) diatomic molecule M- 

These giant Rydberg molecules possess an electronic 
structure characterized by oscillating Born-Oppenheimer 
potential curves evolving from the Rb(n, I > 3) Rydberg 
manifold with binding energies of a few GHz. Their elec¬ 
tronic structure could be controlled and manipulated by 
weak electric fields of a few V/cm [^. The Rydberg- 
field-induced coupling gives rise to a strong hybridization 
of the angular motion of the diatomic molecule together 
with a strong orientation and alignment. In particular, 
the orientation of the diatomic molecule within the Ryd¬ 
berg molecule changes sign and strongly depends on the 
dominant contribution to the Rydberg electric field Qi- 
The two internal rotational states of opposite orienta¬ 
tion of the polar diatomic molecule could be coupled by 
means of a Raman scheme [3] , which might allow for the 
switchable dipole-dipole interaction needed to implement 
molecular qubits 0 . 

Recently, we investigated electronic properties of the 
polyatomic Rydberg molecule composed of a Rb Rydberg 
atom and the KRb diatomic molecule . Our theoretical 
description has included an explicit treatment of the an¬ 
gular degrees of freedom of the diatomic molecule within 


the rigid-rotor approximation. Here, we extend this 
study on the Rb-KRb triatomic Rydberg molecule. We 
investigate the molecule and obtain Born-Oppenheimer 
potential (BOP) curves of the Rb(n, I > 3)-KRb Ryd¬ 
berg molecule with varying principal quantum number 
n of the Rydberg electron and varying internuclear dis¬ 
tance between the ionic core Rb'*' and the diatomic KRb. 
The lowest-lying BOPs evolving from Rb(n, ^ > 3) and 
KRb(A^ = 0) present potential wells with depths of a few 
GHz. We also explore the adiabatic potentials of the Ry¬ 
dberg molecule formed by the KRb being in a rotational 
excited state, and the Rb(26d) Rydberg manifold. These 
BOPs evolving from the Rb(26d) state present wells with 
depths of a few MHz that support several vibrational 
bound states. This opens the possibility of creating these 
macroscopic Rydberg molecules by two-photon excita¬ 
tion of ground-state Rb in an ultracold mixture of Rb 
and KRb. The orientation and alignment of the diatomic 
molecule within the Rydberg molecule is also analyzed in 
terms of the contributions to the electric field due to the 
Rydberg electron and ionic core. 

The paper is organized as follows. In Isection HI we 
describe the adiabatic Hamiltonian of the triatomic Ry¬ 
dberg molecule. The electronic structure of a Rb-KRb 
Rydberg molecule is analyzed in detail in Isection Hll 
We analyze the Born-Oppenheimer potentials evolving 
from the degenerate Rydberg manifolds Rb(n, / > 3) 
and Rb(26d), and the directional properties of the KRb 
diatomic molecule within the Rb-KRb Rydberg trimer. 
The conclusions are provided in Isection IVI 

II. THE ADIABATIC MOLECULAR 
HAMILTONIAN 

We consider a triatomic Rydberg molecule formed by 
a Rydberg atom and a diatomic heteronuclear molecule, 
which is located on the Z axis in the laboratory fixed 
frame (LFF) at a distance R from the ionic core, the 
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latter being placed at the origin of the LFF. The Born- 
Oppenheimer Hamiltonian reads 

Had = Ha+ Hmol ( 1 ) 

where Ha is the single electron Hamiltonian describing 
the Rydberg atom 

with Vi{r) being the /-dependent model potential [H, 
and I the angular quantum number of the Rydberg elec¬ 
tron. 

The term H^oi is the Hamiltonian of the rotational 
motion of the diatomic molecule in the electric field due 
to the Rydberg electron and ionic core 

iJw = HN2-d-F,yd(R,r) (3) 

with N being the molecular angular momentum operator, 
B the rotational constant, and d the permanent electric 
dipole moment. Note that we describe the rotational 
motion of the diatomic molecule within the rigid-rotor 
approximation. The electric field created by the ionic 
core and Rydberg electron at position R reads Hi 

Fryd{^, r) = ® |j. _ J ^|3 

with R = RZ. We consider all three spatial components 
of this electric field taking into account that the position 
of the polar molecule is fixed along the LFF Z-axis. For 
a detailed expression of this electric field, we refer the 
reader to Refs. BP- 

The orbital angular momentum of the Rydberg elec¬ 
tron and the rotational angular momentum of the di¬ 
atomic molecule are coupled to the total angular mo¬ 
mentum of the Rydberg molecule (excluding an overall 
rotation) J = N -|- 1. For the considered configuration 
of the Rydberg molecule, the eigenstates are character¬ 
ized by the magnetic quantum number Mj, which is the 
projection of J onto the LFF Z-axis. 

The Schrodinger equation associated with the Hamil¬ 
tonian [T] is solved by a basis set expansion in the coupled 
basis 

mi—I Mn—N 

fpnlNJMj (r, ^d) = E E {ImiNMNlJMj) 

mi —— I Mn ——N 

^'4^nlm(j^'} ^ ^NMn (5) 

where ipnimij) is the Rydberg electron wave function, 
YnMn i^d) the spherical harmonics representing the field- 
free rotational wave function of the diatomic molecule 
with the rotational and magnetic quantum numbers N 
and Mtv, respectively, fid the internal coordinates of the 
diatomic molecule, and {lmiNM]y\JMj) the Clebsch- 
Gordan coefficient. Taking into account the azimuthal 


symmetry of the Rydberg molecule, the basis set expan¬ 
sion of the wave function reads 

'I'(r, Hd; i?) = ^2 CnlNjiR)tpnlNJMjiT^,^d), (6) 

n,l,N,J 

where we have explicitly indicated the dependence of 
the wave function and expansion coefficients CniNj{R) 
on the internuclear distance R between the diatomic 
molecule and the ionic core. For computational reasons, 
of course, we have to cut the infinite series (l6|) to a finite 
one. In this work, we have included the rotational exci¬ 
tations with iV < 6 for KRb, and for Rb, the (n, I > 3) 
degenerate manifold, and the energetically neighboring 
levels (n -I- Ijd, (n -|- 2)p, and (n -|- 3)s. The quantum 
defect of the Rb(n/) Rydberg state has been neglected. 

III. THE ELECTRONIC STRUCTURE OF THE 
RB-KRB RYDBERG MOLECULE 

In this section, we investigate the electronic struc¬ 
ture of the Rydberg molecule formed by a Rb atom 
and the heteronuclear molecule KRb with rotational con¬ 
stant B = 1.II4 GHz [T^, and electric dipole moment 
d = 0.566 D [l^. The interaction between the electric 
field created by the ionic core and Rydberg electron and 
the permanent electric dipole moment of KRb is respon¬ 
sible for the binding mechanism of Rb-KRb. The electric 
field (U) decreases as the distance between Rb+ and KRb 
is increased, and, therefore, the coupling between the 
subsystems Rb and KRb also decreases. For large enough 
values of R, the adiabatic electronic potentials approach 
the energies of the uncoupled system Eni + BN{N + I), 
with Eni being the energy of Rb(n, /), and BN{N+1) the 
rotational energy of KRb in an eigenstate with rotational 
quantum number N. 


A. The adiabatic potentials evolving from the 
Rb(n, I > 3) degenerate manifold 

The BOPs of the Rb-KRb molecule evolving from the 
Rydberg degenerate manifold Rb(n = 25, / > 3) and with 
magnetic quantum numbers Mj = 0 and Mj = I are 
presented in [Figure l| (a) and [Figure ij (b), respectively. 
In this region of the spectrum, we encounter adiabatic 
states evolving from Rb(28s) and KRb in a rotationally 
excited state, which show a very weak dependence on 
the internuclear distance. This behavior is due to the 
small impact of the Rydberg electric field on the N = 6 
excited state of KRb, which possesses a large field-free 
rotational energy. The electronic states evolving from 
the Rb(n = 25 ,1 > 3) manifold are shifted in energy 
from the uncoupled-system energy Eni + BN{N + I). 
These energy shifts strongly depend on the internuclear 
separation of Rb-KRb due to the interplay between the 
electric fields created by the Rydberg electron and the 
ionic core. They are of the range from a few to tens of 
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FIG. 1. For the Rb-KRb Rydberg molecule, the Born-Oppenheimer potentials as a function of the distance between KRb and 
Rb+ for (a) Mj = 0 and (b) Mj = 1 are shown, specifically the BOP evolving from Rb(28s) and KRb in a rotational excited 
state, and from the Rb(n = 25 ,1 > 3) manifold and the rotational ground state of Rb. The energetically lowest-lying adiabatic 
potentials evolving from the Rb(n = 25, Z > 3) manifold are shown with thicker lines. The zero energy has been set to the 
energy of the Rb(n = 25 ,1 > 3) degenerate manifold and the KRb {N = 0) state. 


GHz. The BOPs oscillate as a function of R, and nu¬ 
merous avoided crossings among neighboring electronic 
states characterize the spectrum. For R ~ 1000 ap, the 
energetically lowest-lying BOPs present potential wells 
of a few GHz depths, which support several vibrational 
states 0. By further increasing R, the coupling between 
the Rydberg atom and KRb decreases, and the BOPs 
approach the asymptotic limit of the uncoupled system. 
Note that in these figures, the zero energy has been set 
at i?n= 25 ,i> 3 , i- e., for the uncoupled system KRb(N = 0) 
and Rb(n = 25,1 > 3). 

We focus now on the evolution of the electronic struc¬ 
ture as the principal quantum number of the Rydberg de¬ 
generate manifold Rb(n, Z > 3) is increased. In [Figure 2[ 
we present the energetically lowest-lying adiabatic elec¬ 
tronic states with Mj = 0 and Mj = 1 evolving from 
the Rb(n, I > 3) manifolds for n = 25,26, 27, 29, and 31. 
These BOPs show a qualitatively similar oscillatory be¬ 
havior as a function of R. By increasing the principal 
quantum number of the Rydberg manifold, some differ¬ 
ences are observed. First, the energy shift of these BOPs 
from the degenerate manifold of the Rydberg atom de¬ 
creases as n is increased. Second, the location of the 
potential wells is shifted towards larger values of the in- 
ternuclear separation between Rb'*" and KRb. This is due 
to the radial extension of the electronic Rydberg wave 
function of Rb(n, I > 3) which increases as n is enhanced. 
Third, the depth of the wells decreases with increasing n. 
The outermost potential well of the Mj = 0 BOPs is very 
shallow and do not support vibrational bound states. In 
contrast, the potential well with Mj = 0 located at the 
left of this outermost one is the deepest one supporting 
several bound states and with depths of several GHz for 
all the Rb(n, I > 3)-KRb(N = 0) molecules investigated 
here. For instance, we have found 7 and 6 vibrational 
bound states in these potential wells of these energetically 
lowest-lying Rb(n = 25,1 > 3) and Rb(n = 31, Z > 3) 


BOPs with Mj = 0, respectively. 

We analyze now the directional features of the 
KRb molecule within the Rb-KRb triatomic Rydberg 
molecule. For the lowest-lying electronic potentials 
evolving from the Rb(n = 25 ,1 > 3) and Rb('n = 29 ,1 > 
3) manifolds and with Mj = 0, the orientation, (cosOj), 
and alignment, (cos^dd), of KRb along the LFF Z-axis 
are shown in |Figure~3] (a) and [Figure 3| (b), respectively. 
For the sake of simplicity, we only analyze the results 
for the Rb-KRb in these two manifolds, but qualitatively 
similar behavior is obtained for the others BOPs pre¬ 
sented in [Figure 2] We have computed the orientation 
and alignment of the rotational ground state of KRb in 
an electric field of varying held strength 5.14/i?^ GV/cm, 
i. e., the strength of the electric held due to the Ry¬ 
dberg core, being parallel and antiparallel to the LFF 
Z-axis, see [Figure 2| These results allow us to get a 
deeper physical insight into the interplay between the 
electric held due to the Rydberg electron and the Rb+ 
core. For both BOPs, we encounter a regime where the 
electric held due to the Rydberg core is dominant, i. e., 
the total electric held is antiparallel to the LFF Z-axis, 
and KRb is anti-oriented. In this regime, the orienta¬ 
tion of KRb within Rb-KRb is very close to the orienta¬ 
tion of KRb in the 5.14/i?^ GV/cm electric held. Let 
us comment, that for R = 1000 oq, the electric held 
strength due to the Rydberg core is 5 kV/cm. The os¬ 
cillations of the orientation of the KRb within Rb-KRb 
are due to the Rydberg electrons electric held, which 
oscillates resembling the behavior of its wave function. 
The amplitude of these oscillations increases as the Ry¬ 
dberg electrons electric held becomes more important. 
When the contribution of this held becomes dominant, 
the orientation of KRb changes sign and KRb becomes 
oriented along the Z-axis. The differences in the elec¬ 
tronic Rydberg wave function of Rb(n = 29 ,1 > 3) and 
Rb(n = 25, Z > 3), are rehected in the fact that the KRb 
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FIG. 2. For the Rb-KRb molecule, the energetically lowest-lying states with Mj = 0 and Mj — 1 evolving from the Rb(n, I > 3) 
Rydberg manifold, for n = 25, 26, 27, 29, and 31, and KRb in the rotational ground state N = Q are shown. To facilitate the 
comparison, the zero energy has been set to the energy of the Rb(n, I > 3) degenerate manifold and the KRb(77 = 0) level. 
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FIG. 3. For the KRb molecule within the energetically lowest-lying BOPs with Mj = 0 of Rb(n = 25 ,1 > 3)-KRb (thin solid) 
and Rb(n = 29 ,1 > 3)-KRb (thick solid), we show (a) the orientation {cos Od) and (b) the alignment (cos^ 9d) along the LFF 
Z-axis versus the internuclear distance between Rb’*' and KRb. We also present the orientation and alignment of KRb in an 
electric field parallel (thick dotted) and antiparallel (thin dotted) to the LFF Z-axis and strength 5.14/7?^ GV/cm. 


molecule is anti-oriented for a larger range of values of R 
for Rb(n = 29, ^ > 3)-KRb(A^ = 0). KRb is maximally 
oriented around the classical turning point R 2n^ao, 
and we have {cos 9d) ~ 0.8, which is the orientation of 
KRb in a strong electric field of 50 kV/cm strength M- 
The internuclear separation of Rb-KRb at which the KRb 
orientation reaches the maximum is shown versus the 
principal quantum number of the Rb(n, I > 3) degener¬ 
ate manifold in |Figure 4| We observe a very good agree¬ 
ment with the 2n^ao behavior. In [Figure 4] we show also 
the internuclear separatioir of Rb-KRb at which the en¬ 
ergetically lowest-lying BOP with Mj = 0 reaches the 
two outermost minima. The radial position of these two 
minima is shifted to lower values of R compared to the 
maximum of the KRb orientation. By further increasing 
R, the orientation decreases towards zero and approaches 
the orientation of KRb in the 5.14/i?^ GV/cm electric 
field. Note that in the absence of an electric field, the 
KRb molecule is not oriented. 

The alignment along the LFF Z-axis of KRb within 
Rb-KRb also shows an oscillatory behavior around the 
alignment of the KRb in a 5.14/i?^ GV/cm electric field. 
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FIG. 4. For the Rb-KRb Rydberg molecule, we show the 
distance between the Rb"*" core and KRb (•) at which the 
orientation of KRb is maximal within the lowest-lying Mj — 0 
adiabatic potentials evolving from the Rb(n, I > 3) manifolds 
for varying principal quantum number n. For these lowest- 
lying BOPs, we also present the internuclear separation R of 
Rb-KRb at which the one before the outermost (□) and the 
outermost (■) potential wells reach the minimum value. For 
comparison, we show the classical turning point R ~ 2n^ao 
(*). 
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cf. [Figure 3] (b). As the coupling of the dipole moment 
of KRb with the electric field due to the Rydberg elec¬ 
tron becomes dominant compared to the coupling with 
the Rb+ electric field, the amplitude of the oscillations 
become very large. When the distance of KRb form the 
Rb’*' core is large, the alignment approaches the field-free 
value (cos^ 6d) ~ 0.3333. 

B. The adiabatic potentials evolving from the 
Rb(26d) Rydberg state 

In the absence of an external electric field, the 
ultralong-range Rydberg molecule could be created from 
an ultracold mixture of Rb and KRb, by exciting the Rb 
from the ground-state to a nd Rydberg level in a two- 
photon process. In this section, we investigate the BOPs 
of the Rydberg molecule formed from Rb(26d). 

The adiabatic potentials with Mj = 0 evolving from 
the Rb(26d) Rydberg states and the diatomic molecule 
in the rotational excitations N = 2, 3 and 4 are present 
in [Figure 5| There are 5 electronic states with Mj = 0 
and total angular momentum J = 2 — N^. N -\- 2 in 
the uncoupled system. These BOPs show either a strong 
or a weak dependence on i?, and when R is increased 
beyond a certain value they approach the energetical 
limit AE 2 ed + BN{N -I- I), with N = 2, 3 and 4 and 
AE 2 Gd = 1^24,3 — E 2 Gd\- The BOPs with oscillatory be¬ 
havior possess potential wells having depths of a few MHz 
or tens of MHz. These potential wells are deep enough 
to accommodate at least one vibrational level in which 
the polyatomic Rydberg molecule would exist. For a few 
vibrational levels, the square roots of the wave functions 
are shown in the potential wells in [Figure 5] For the 
BOPs evolving from the N = 3 and = 4 rotational 
states of KRb, we encounter several vibrational bound 
states with wave functions extending over several poten¬ 
tial wells. 

In [Figure 5] (d), we present the orientation of the di¬ 
atomic molecule along the LFF Z-axis for the electronic 
potentials evolving from the Rb(26(i) Rydberg manifold 
and KRb {N = 3). The diatomic molecule is weakly ori¬ 
ented along the LFF Z-axis, and similar results are ob¬ 
tained for the BOPs evolving from the diatomic molecule 
in the rotational excitations N = 2 and N = A. This 
weak orientation can be explained in terms of the small 
impact of the Rydberg electric field on the rotational 
excited state, which is smaller than for the rotational 
ground state due to its larger rotational energy M- 
Thus, the weaker the coupling between the N rotational 
state and the neighboring levels iV ± 1, the smaller the 
orientation is. 


IV. CONCLUSIONS 

In this work, we have considered the ultralong-range 
Rydberg molecule formed by the KRb diatomic molecule 


and the Rydberg Rb atom. The coupling mechanism of 
this triatomic Rydberg molecule is based on the interac¬ 
tion between the electric field produced by the Rydberg 
electron and core and the permanent electric dipole mo¬ 
ment of the polar molecule. We have performed an anal¬ 
ysis of the electronic structure of Rb-KRb when the Rb 
atom is in the Rydberg degenerate manifold {n,l > 3), 
with n = 25,26,27,29, and 31. These adiabatic poten¬ 
tials show an oscillatory behavior as a function of the 
separation between Rb+ and KRb. For the Rb(n, / > 3)- 
KRb(V = 0) molecules, the lowest-lying BOPs show po¬ 
tential wells of a few GHz depths. As the principal quan¬ 
tum number of the Rydberg excitation n is enhanced, the 
depths of these potential wells decrease and their posi¬ 
tions are shifted towards larger values of the internuclear 
separation R. These features reflect the oscillatory struc¬ 
ture of the Rydberg electronic wave function. 

We have also investigated the properties of the Ry¬ 
dberg molecule formed by the Rb(26(i) Rydberg state 
and the diatomic molecule in an excited rotational state. 
BOPs with both an oscillatory and a smooth behav¬ 
ior as function of R are encountered. For Rb(26(i)- 
KRb(N), these oscillatory BOPs present potential wells 
with depths of a few tens of MHz, which can accom¬ 
modate a number of bound vibrational states. These 
results demonstrate that these ultralong-range Rydberg 
molecules could be created by a standard two-photon ex¬ 
citation of ground state Rb in an ultracold mixture of Rb 
and KRb. 

Due to the electric field of the Rydberg electron and 
core, the directional properties of the KRb molecule 
within the Rb-KRb triatomic molecule are significantly 
affected. The orientation and alignment of KRb has been 
analyzed in terms of the contribution of the electric fields 
due to the Rydberg electron and ionic core. Within these 
BOPs, the polar diatomic molecule presents two opposite 
orientations along the LFF Z-axis with its electric dipole 
moment pointing either toward or away from the ionic 
core. 
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FIG. 5. For the Rb-KRb Rydberg molecule, Born-Oppenlieimer potentials with Mj = 0 evolving from the Rb(26d) state and 
KRb in the rotational levels (a) N = 2, (b) N = 3, and (c) N = 4. We have plotted the 5 electronic states with Mj — 0 and 
total angular momentum J = 2 — N,...N-\-2,m the uncoupled system. In the potentials wells, the square of the vibrational 
wave functions are shown in arbitrary units. The zero energy has been set to the energy of the Rb(n = 25 ,1 > 3) degenerate 
manifold and the KRb {N = 0) state. In panel (d), we show the orientation along the LFF Z-axis of KRb within the BOPs of 
the Rydberg molecule Rb(26d)-KRb(A^ = 3). 
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